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ABSTRACT 

We examine the vertical structure and equilibrium of the molecular gas layer 
in the galactic disk, measuring its scale height and velocity dispersion as a 
function of Galactic radius by modeling the CO emission at the tangent points. 

The model takes into account emission from a large path length along the 
line of sight, corresponding to an interval (AR) of typically a few hundred 
parsecs in galactic radius; and is parametrized by the scale height of the gas, the 
centroid in z, the rotation velocity and the velocity dispersion; these parameters 
are assumed to be constant over the interval AR. This model is then fit to the 
12 CO survey of Knapp et al.(1985) to determine the best fit parameters. 

The terminal velocity values are found to be in good agreement with those 
obtained from HI data. The x 2 P er degree of freedom from fitting the models 
(considering only the photon statistics) range from 2 to 19, with a median value 
of 6. The main source of error is found to be the 'shot noise' due to the small 
number of clouds. Simulations of the tangent point emission using discrete 
molecular clouds are carried out to estimate the errors. We use the observed 
radial distribution of molecular gas and the 'standard' size-linewidth relation for 
molecular clouds in the simulations. Modeling the simulations gives values of x 2 
per degree of freedom similar to those obtained from modeling the observations. 

The scale height of the distribution is found to increase fairly monotonically 
with radius. The value of the velocity dispersion varies from ~ 2kms _1 to 
llkms^ 1 with a typical uncertainty of ~ 3 kins" 1 . The variation in velocity 
dispersion is consistent with a monotonic increase with galactic radius. 

The midplane mass density of the disk po{R) is calculated from the 
scale-height and velocity dispersion (under the assumption that the velocity 
dispersion is isotropic) and is consistent with the local value po(Ro) of 
O.2M pc -3 determined from stellar kinematics. 

Subject headings: interstellar gas, Galaxy, internal motions, mass 
models 
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1. INTRODUCTION 

Because of its expected simplicity, the distribution and kinematics of gas in the 
direction perpendicular to the plane of the disk of the galaxy should have much to tell us 
about the forces acting on it and its state of equilibrium. Considering the most obvious 
force, gravity, and recalling that the scale height of molecular gas is much smaller than that 
of stars, we expect a Gaussian distribution in z p(z) = p exp(—z 2 /2a 2 ) (Spitzer 1942), 
if the velocity dispersion a v of the gas is constant. The Gaussian scale height is then 
a z = °"u/4vr(jpo) where po is the mid-plane mass density of the disk. 

Knowing the scale height and the velocity dispersion in the vertical direction one can 
estimate po, the mid-plane mass density of the disk. Studies using stellar populations (F 
dwarfs and K giants) as test particles to probe the potential give estimates of both the local 
disk mass density po an d the integral surface mass density S (Bahcall et al. 1992, Bahcall 
1984a, Bahcall 1984b, Bahcall 1984c, Bienayme et al. 1987, Kuijken & Gilmore 1989a, 
Kuijken & Gilmore 1989b, Kuijken & Gilmore 1989c, Kuijken & Gilmore 1991a, Kuijken 
& Gilmore 1991b, Kuijken 1991). In a complimentary approach we can hope to probe the 
potential at different galactocentric distances with gas (clouds) since gas emission from 
the whole disk can be observed. Since the gas disk is much thinner than the stellar disk, 
its vertical distribution probes only the midplane mass density. A similar study has been 
carried out for HI by Celnik et. al (1979) and Merrifield (1992). 

The kinematic information about molecular clouds, especially their velocity dispersion 
is important for determining their dynamics, with each other as in collisions (influencing 
cloud formation/destruction, star formation); and with the stars, increasing the velocity 
dispersions of stars (Spitzer & Schwarzschild 1951). The variation of velocity dispersion 
with Galactic radius may give clues to the sources of random kinetic energy to the clouds 
(Gammie et al. 1991), star/cloud formation as a function of galactic radius etc. 

Since the molecular gas forms a cold, thin layer in the disk, it may be possible to look 
for signatures of disequilibrium in departures from Gaussianity of the distribution. These 
departures may be due to violent events stirring up the gas (for example supernovae), short 
lifetimes of molecular clouds, infalling gas etc. or because of more than one isothermal 
population with different values of velocity dispersion a v . 

The large scale distribution of molecular gas has been mapped using the Massachusetts- 
Stony Brook galactic plane survey of 12 CO in the plane by Sanders et al. 1986, Clemens et 
al. 1986; and the Columbia survey by Cohen et al. 1986, Bronfman et al. 1989, Dame et 
al. 1987. We use data from the Bell Labs 12 CO survey (Knapp et al. 1985), consisting of 
strip maps in latitude because we need good latitude coverage and sacrifice close sampling 
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in longitude. The details of the survey and the data as they relate to the present analysis 
are presented in section 2.1. 

The scale height of the molecular gas has been previously estimated by Clemens et al. 
(1988) to be about 65 pc at R = 0.95i? and by Sanders et al. (1984), Bronfman (1988), to 
be about 60 pc. Different surveys thus give very consistent values of the scale height. 

The same cannot be said for the velocity dispersion a v of the clouds. The values of 
a v found by different groups are inconsistent by a factor of two, if not more. The velocity 
dispersion for local (distance from the sun < 3kpc ) clouds has been estimated with the 
additional input of their distances and the Galactic rotation curve by Stark (1984) and 
Stark and Brand (1989), to be 7.8lg;| kms' 1 . Blitz et al. (1984) find a v = 5.7 ± 1.2 kms" 1 
for high latitude (hence local) clouds. Liszt and Burton (1983) estimate a single value of 
a v = 4.2 km s" 1 for the inner galaxy by measuring the dispersion of what they call 'jitter', 
which is the difference between the terminal velocity of HI and that of 13 CO. Clemens 
(1985) estimates a similar value by fitting a Gaussian to the tangent point emission. 
To get rid of emission coming from galactic radii smaller than the tangent point radius 
(sub-tangent point emission) the fit is restricted to ±a z since the emission from nearby 
gas extends to higher latitudes. The present paper takes a new approach to modeling the 
tangent point emission by fitting simultaneously the terminal velocity, velocity dispersion 
and the scale height to the observations. We do not try to isolate the tangent point emission 
but take into account emission from the sub-tangent point gas. The modeling procedure is 
described in sections 2.2 and 2.3. 

We use this model to get the best fit values of the aforementioned parameters for each 
longitude observed. As expected the main source of uncertainty in the parameters is due 
to the dumpiness of the molecular gas distribution. To estimate these uncertainties in the 
parameters measured, we carry out Monte Carlo simulations of the cloud distribution in 
the galaxy, using the known radial distribution of gas and standard cloud properties. The 
details of the simulations are discussed in section 2.5. 

In section 3 we present the best fit parameters, and try to interpret their variation 
across the inner galaxy. The errors and systematic shifts in these parameters are estimated 
from the Monte Carlo simulations. We also estimate the midplane mass density of the disk. 

Section 4 contains a discussion of different values, methods and the definitions of the 
parameters characterizing the galactic distribution as estimated by previous studies; and 
how consistent various determinations are, particularly with respect to velocity dispersion 
estimates. The assumptions behind this study and various caveats are also touched upon. 
The conclusions of this paper are summarized in section 5. 
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Appendix A contains contour maps of CO emission as a function of latitude and 
velocity, for several different longitudes. Contours of emission according to the best fitting 
models are superposed. 

The value of the circular velocity at the Sun 6 is taken to be 220 km s -1 , and the 
assumed distance to the galactic center is 8.5 kpc. 

2. METHOD 
2.1. Data 

The observations analyzed here were made using the 7 meter antenna at Bell 
Laboratories. The data and details of the survey and the instruments are presented by 
Knapp et al. (1985). Here we summarize the information necessary for the present analysis. 

The survey consists of 38 strip maps (in latitude) of the 12 CO(l — > 0) line emission. 
The maps are sampled at intervals of Ab = 2' with a half-power beamwidth of 100". The 
latitude coverage is ~ ±2°, varying slightly from one line of sight to another. The latitude 
extent is more than adequate for studying the tangent point emission (except possibly for 
I > 60°), extending more than three scale heights both above and below the centers of 
distribution in more than half the lines of sight. 

The strip maps were taken at longitudes between 4° and 90°, spaced at equal intervals 
of A sin / = 0.025. The velocity resolution is 0.65 km s -1 . The rms noise is typically 0.3 K 
but varies slightly in each latitude strip. Linear baselines were removed from individual 
spectra. Exceptions to this treatment are few and are listed by Knapp et al. (1985). 

We are able to analyze 22 of the 38 observed lines of sight. The low longitude 
(4° < / < 17°) maps were unusable due to the lack of emission at (or reasonably near) 
the expected tangent point velocities. The high longitude (I > 61°) maps were unusable 
because their tangent point emission is at low velocities and is contaminated by emission 
in reference positions; and because the latitude extent of the maps is not adequate for this 
(almost) local emission. 

2.2. Tangent points 

Assuming a circularly symmetric model of the galaxy we can obtain distances to the 
emission at extreme velocities in the first and fourth quadrants. At each galactic longitude 
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/, extreme velocity emission comes from the tangent point at a galactocentric distance of 
R = R sinl (Figure 1). The tangent points are at a distance d = R cosl from the sun. 
Thus the observed angular extent b gives the height z — Rq cos / tan b. In the first quadrant 
the velocity of the gas at radius R is given by 

VW = !^-e sinZ (1) 

where 6 is the circular velocity at R. At the tangent points equation (1) reduces to 
V T = 6 - O sin/. 

From the observed b extent and the terminal (extreme) velocities one could calculate 
the scale height and the rotation velocity. The velocity profiles however do not have a sharp 
cutoff due to the velocity dispersion of the gas. Ideally the emission from the tangent point 
T(b,v) is a bivariate Gaussian in altitude and velocity; but emission from nearby radii is 
not well separated in velocity because of the velocity dispersion of the gas. 



2.3. The velocity dispersion 

The finite velocity dispersion has two main effects that the present analysis needs to 
take into account. First it makes the velocity-to-distance conversion fuzzy, so it is difficult 
to separate emission from R (at the tangent point) and that from R'(R' — R + AR), where 
AR depends on the velocity dispersion. For cold gas (zero velocity dispersion) a difference 
in velocity AV T = 7kms _1 would correspond to AR ~ 270/ sin(/) pc near the tangent 
point. Expecting velocity dispersions of this size, it is meaningful to calculate parameters 
like the scale height, rotation velocity and the velocity dispersion, averaged over intervals 
AR of a few hundred parsecs. 

At the same time we must take into account the fact that the emission at velocity V 
comes from R, and from the annulus between R and (R + AR), so that T(b,v) is no longer 
a simple bivariate Gaussian. Celnik et al.(1979) derived the expression for the expected 
line shape near terminal velocity. In this section we calculate the 2-dimensional (in latitude 
and velocity) profile of emission from near the tangent point, i.e. between R and R + AR, 
taking the scale height, rotation velocity and velocity dispersion to be constant over AR. 
We expect the additional information in the apparent latitude-extent as a function of 
velocity to yield more constraints on the positions of the emitting regions. 

Consider emission seen at velocities close to the terminal velocity; (AV <C 0) 



V(R') = V T -AV = 9- o sinZ - AV 



(2) 
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giving 



R' Gsin/ 



(3) 



R B-AV 

An observer's line of sight intersects the annulus at radius R 1 at two points (figure 1). The 
distances to the points of intersection, r\ and r 2 are given by 
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For the same scale height, gas at the subtangent points will have smaller and larger 
latitude extent than the tangent point gas, corresponding to the far and the near part of 
the annulus at R' 

(7.1 (V) = C ° Sj , V,(V T ) 

2 cos/ =F smlyJ2AV/G 

Finally we must take into account velocity crowding effects. The optical depth per 
velocity interval is determined by the line of sight distance per unit velocity interval. 
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Thus the emission in the velocity interval AV (corresponding to AR = Q /R sin(l)) 
near the terminal velocity, keeping parameters, a v , a z fixed over the velocity interval is 



t(v,z)=ay: 
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(7) 



2.4. Fitting 



The model profile is calculated (eqn. 7) as a function of latitude b and velocity V 
for each longitude and a least-squares fit is made to the data to determine the best fit 
parameters, namely the amplitude A, the centroid in the vertical direction zq, the vertical 
scale height a z , the terminal velocity Vp and the velocity dispersion a v . 
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The fitting is done by least squares minimization using the downhill simplex routine, 
'amoeba' (Press et al. 1993). Contour maps of emission as a function of latitude b and 
velocity V and superposed best-fit models are given in Appendix A (Figure 12) and Figures 
2 and 3a . 

The velocity range over which the fit is made is determined by the width of the extreme 
velocity feature. For each line of sight the spectra at all latitudes are summed to form a 
composite spectrum. The peak at the highest velocity is identified as the terminal velocity 
feature. The lower velocity at which the emission drops to half- maximum is defined as Vhaif- 
We fit the emission seen at velocities greater than Vhaif. 

To see how sensitive the parameters obtained are to a change in the velocity range 
over which the fitting is done, we do the fitting for different values of Vhaif for one of the 
longitudes / = 39. Figure 3a shows the various values of Vhaif in relation to the terminal 
velocity feature. The best fit parameters for each of the fits is plotted against Vhaif m figure 
3b. We see that the parameters are fairly insensitive to the velocity range over which the 
model is fit. 



The above analysis holds true for a diffuse medium. For molecular clouds we might 
interpret T(b, V) as a smoothed probability of finding molecular clouds at galactic latitude 
b and velocity V. The molecular gas is distributed as clouds and cloud complexes, which 
sample this probability. The main uncertainty in deriving the parameters of the models is 
the shot noise due to the small number of clouds and the dumpiness of the clouds. From 
figures 2, 3 and 12, it is clear that the models do not trace the detailed structure of the gas. 

If we ignore the dumpiness of the molecular gas layer and treat the emission as if it 
were from a smoothly distributed medium, the \ 2 P er degree of freedom is typically found 
to be a few (2-19). The x 2 per degree of freedom, or reduced x 2 is calculated by the relation 



where T v h oton is the photon noise, and N the number of degrees of freedom. Since the main 
source of noise is the dumpiness of the medium, which is greater than the instrumental 
noise, we do not expect the x 2 to be 1. On the other hand we use the relation x 2 = 1 to 
estimate 'realistic' noise 'T rea /' of the emission and then proceed to calculate the confidence 
intervals on the parameters using the rescaled noise. 



2.5. 



Error estimation 




(8) 
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While setting the x 2 — 1> we nave lost information on the goodness of fit of the models. 
Also l T rea i is not necessarily normally distributed, an assumption made while estimating 
confidence levels on the distribution parameters. So we carry out Monte-Carlo simulations 
of the cloud population in the galaxy. This also useful for detecting systematic errors, to 
the extent that we understand and can simulate the cloud population. 

If we assume the clouds to be ballistic particles, it is the locations and velocities of the 
clouds that follow the distribution derived in the last section. Clearly a large number of 
optically thin entities would mimic a smooth distribution. The velocity crowding (dr/dv 
— > oo) at the tangent points somewhat justifies the assumption that there is a large number 
of clouds. Velocity intervals of few kms -1 near the terminal velocity correspond to a 
kpc or so of distance r traversed along the line-of-sight. We would expect ~ 40 clouds in 
the velocity interval of 7kms _1 near the terminal velocity (at longitude I = 35°). As for 
'transparency', molecular clouds show optically thick clumps and filaments with a typical 
surface filling fraction of 30% (e.g. Bally et al. 1987). Thus the approximation that we are 
counting the number of clouds/clumps, and that little emission from clouds is shadowed, is 
perhaps not entirely unjustified. 

To test the limitations of this approximation, and to get an estimate of the confidence 
levels on the parameters, we simulate the cloud population in the Galaxy, taking the total 
mass of the molecular gas in the galaxy to be 2 x 1O 9 M , and the proportionality constant 
X = M§± to be 3 x 10 20 cm- 2 K- 1 . 

With the additional input of the size-linewidth-mass distribution, 

M \ ( Av i n \ 

(9) 




V1OOM ; V 0.36 km s 
and the mass spectrum of the clouds (Solomon & Rivolo 1989) 

dn(M) _ 



dM 



oc M 5 (10) 



we can perform Monte-Carlo simulations of the cloud distribution in the galaxy. The upper 
and lower cutoffs for the cloud mass are taken to be 10 6 M & and 10 3 M & respectively. The 
cloud centers are distributed in (R,9,z). The radial distribution of the gas is taken from 
earlier galactic surveys, which closely sample in longitude (Sanders et al. 1986, Dame et 
al. 1987). The azimuthal distribution at a particular radius is taken to be uniform. In 
the z-direction the clouds centers follow a Gaussian distribution with mean z = Opc, and 
scale height a z = 50 pc. The cloud-cloud velocity dispersion is taken to be a v = 7kms~ 1 . 
The clouds themselves are represented as Gaussians in z and v, with the aforementioned 
size-line width relationship between a z (cloud) and a v (cloud). 



9 



The clouds are then 'observed' with a beam of FWHM 100" for 9 longitudes spanning 
the latitude range of the survey. The simulated data sets have the same latitude and velocity 
coverage as the survey and were modeled in the same fashion. Cloud-cloud shadowing near 
the terminal velocity is considerable (~ 7 for mid-longitudes) and is handled by scaling 
down the contribution from the shadowed region of a cloud by / = (1 — 0.3) n where n is the 
number of clouds between that cloud and the sun. Detailed modeling of the beam shape, 
cloud sizes, radiative transfer etc. has not been carried out. A more complicated model is 
perhaps not justified in view of the severe blending of the features at the tangent point. 

3. RESULTS 

3.1. Results of the simulations 

Examples of contour maps (in b-v) derived from the Monte-Carlo simulations are 
shown in Figure 4 for longitudes I = 22.02°, I = 50.8° and I = 35.1° and the best fit models. 
At 1=35.1° the terminal velocity regions are severely crowded with typically 40 clouds in 
the region we are trying to fit. The cloud-cloud shadowing goes up to 7, i.e. we might be 
looking 'through' up to 7 clouds near the tangent point velocities. 

The parameters resulting from fitting the models to the simulations of 9 longitudes 
are shown in Table II and figure 5. All the longitudes simulated have identical input 
distribution parameters, to test for systematic errors as the velocity crowding changes 
with longitude. The parameters measured agree well with the input parameters, with the 
exception of the scale height a z , which is measured to be higher than the input. It is 
reassuring that the x 2 P er degree of freedom for the simulated models is found to vary 
between 6 and 13, similar to the values for the observations (2-19). The velocity dispersion 
a v has an uncertainty of ±3kms _1 , while the error bars obtained from leastsquare fitting 
and scaled up noise correspond to ±3.5 km s" 1 . The error-bars for a z from Monte-Carlo 
simulations are larger than the ones from the leastsquare analysis. There are no systematic 
trends in the estimates of the parameters with longitude. 

The scale height of the gas is expected to be higher than the scale height of the cloud 
centers, since o" cente rs and a c i ou d add in quadrature to give the scale height of the gas. We 
repeat the simulations with clouds of different sizes at 1=35.1° , and the measured scale 
height of the gas does go down with the size of the clouds. The values of a z however 
are still slightly higher than expected if the 0" centers and <7 c i ou d add in quadrature. This 
is because of crowding/shadowing which tends to flatten the z-distribution. Simulations 



10 



with no shadowing (i.e. transparent clouds) and all of the mass in 10 4 M and 10 3 M Q 
clouds reproduce the expected scale height of the gas (cf Table III). For the standard 
mass-spectrum of the clouds in 10 6 M to 10 3 M & clouds the contribution to scale height 
defined by a^(cloud) = ^(measured) — of(centers) = 56 pc. Only one longitude in the 
Knapp et al. survey shows a scale height larger than 55 pc (table I), making it difficult to 
accommodate the size spectrum used for the simulations. 

The scale heights estimated for the Monte-Carlo simulations depend on other properties 
like the density profiles of the clouds. Using clouds with a sharp cutoff in emission at the 
cloud boundaries as seen in local clouds ( e.g. Blitz & Thaddeus 1980), the measured 
scale height is lower than estimated for clouds with Gaussian density profiles (table 
II). Simulations of emission at 1=35.1° were done with different values of scale height 
^(centers) = 2, 10, 20, 30 and 40 pc, and a size spectrum of 10 3 — 10 6 M e . The contribution 
of the cloud size to the scale height a^(cloud) is constant ~ 42 ± 2.7 pc (figure 6). The scale 
height of molecular gas at the tangent points is observed to be < 40 pc along many lines of 
sight. If the resulting scale height of the simulated data sets were not sensitive to cloud 
properties we would conclude that this analysis favors small clouds. 

Since the molecular gas is clumped, we wish to investigate the effect of having no gas 
at the tangent points. It has been shown by Shane & Beiger-Smith (1966) that tangent 
points with little or no gas show a low, broad-shouldered profile (also see figure 7). For the 
distribution of clouds a lower profile may be taken to mean that statistically there will be 
fewer clouds at the extreme velocities. Since we allow the fitting program to choose the 
velocity dispersion a v and the terminal velocity Vr, missing tangent point gas will lead to a 
lower value of V T and a higher value of a v . The other possibility is that the fitting routine 
will get confused by the small number of clouds and fit a single cloud giving a small a v . 
The scatter in the parameters is found to be fairly large for the simulations with no gas at 
the tangent points, as expected from a smaller number of clouds (figure 8). Thus the very 
high (llkms -1 ) as well as the very low (2 kms -1 ) values of a v measured from the data are 
not believable. 

3.2. Distribution parameters 

So far we have derived the general functional form of T(b,V) and the shape of the 
emission contours. We also see that the derived parameters are robust to small changes in 
the fitting procedure, and the derived uncertainties in the parameters are comparable to 
those from the simulations. All the data and the best fit models are shown in Appendix A, 
figure 12. 
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Figures 9 and 10, and table 1 show the distribution parameters Vp, cr v , a z , z for all the 
longitudes. At low longitudes (1=17°, 19°) a v and a z are found to have low values but with 
larger error-bars, due to the relatively small amount of CO at those radii. 

The midplane position zq as derived from the fitting shows a smooth undulation (figure 
10). The scale height increases with radius R; also the errors in a z increase, mainly due 
to the coverage in latitude translating to a smaller coverage for the nearby gas at higher 
longitudes. A linear fit gives a gradient in scale height of 5.5 pc kpc -1 . 

The terminal velocity closely matches the terminal velocity curve for HI (Gunn et al. 
1979), pointing to the similarity of the HI and CO distribution and kinematics. This also 
reassures us that the tangent points are not substantially devoid of molecular gas, since HI 
has been shown to be present at tangent points by the absence of low broad-shouldered 
profiles (Shane and Beiger-Smith 1966). 

The values of the velocity dispersion found for different lines of sight range from 
~ 2kms~ 1 to 11 kms -1 . There is however a trend for the a v to increase at higher longitudes, 
i.e. with increase in galactocentric radius. A straight line fit through the a v values shows 
an increase from 3.8 kms -1 at 2.5 kpc to 7.1 kms -1 at 7.5 kpc. 



3.3. Midplane mass density of the disk 



Considering the condition for vertical equilibrium of the clouds in the galactic disk, 
in the simplest case of a single isothermal population of clouds, the vertical component 
of the gravitational force is balanced by the pressure due to turbulent motions of the gas 
Pz = Pz&v- The thickness of the molecular gas layer being much less than that of the mass, 
the vertical component of the force K z = —AirGp^z (po is the midplane mass density) and 
the distribution of gas is a Gaussian with scale height a z . The midplane mass density of the 
disk is then given by 



With the additional assumption of isotropy of velocity dispersion of clouds, we can 
estimate p from the scale height and the velocity dispersion (Knapp 1988). 

Figure 11 shows the midplane mass density of the disk computed from equation 10 
and the best fit parameters oy and o z for different longitudes. We fit an exponential disk 
p = p c exp(— R/R c ), yielding a scale length R c = 5pc and midplane mass density at the 
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center of the disk po(0) = 11 M Q pc~ 3 . From these two values we get the midplane mass 
density at sun po(Ro) — O.2M pc -3 . The scale lengths of 3.2 kpc and 10.1 kpc and central 
densities of po(0) = O.5M pc~ 3 and 0.07M Q pc~ 3 lie within 66% confidence limits, leading 
to an uncertainty of a factor ~ 2.5 in the value of pq(Rq). 

The simulations however show that the parameters determined have systematic errors. 
The scale height of the gas is measured to be systematically higher than the scale height of 
the cloud centers. If we knew the sizes of the molecular clouds reliably we could compensate 
for that. Liszt and Burton (1981) favor a smaller (5-10 pc) typical cloud size by comparing 
the morphology of the 1-v diagrams of simulations with observed data. If most of the mass 
of the molecular is in small clouds the corrections to the scale height are small. Shadowing 
of one cloud by another flattens the distribution and leads to a further overestimate of the 
scale height. This effect should be smaller for higher longitudes (less crowding/shadowing); 
so we probably are overestimating the scale length R c . Apart from these effects seen in 
the simulations, we may expect the velocity dispersion measured to have a non-random 
component from streaming motions as we are measuring the horizontal velocity dispersion. 
Further discussion of this effect is postponed to the next section. 

4. DISCUSSION 

For a consistent analysis of the vertical equilibrium in the disk, we must measure the 
scale heights and the velocity dispersion of the same population. In the current analysis 
we do just that for the molecular gas seen at the tangent points. The tangent point 
emission is modeled in two dimensions, so the scale height and the velocity dispersion are 
fit simultaneously, along with the other two free parameters of the model, the rotation 
velocity and the centroid of the z-distribution. One might worry about correlated errors in 
the determination of the four quantities, for example a lower terminal velocity estimation 
may lead to a higher measurement of the velocity dispersion. While we don't know the real 
parameters we can estimate the local dips and peaks in the parameters as compared to 
the values smoothed over large radii. The fluctuations in different parameters do not seem 
to be correlated (figures 9 and 10). Moreover the best-fit parameters of the Monte-Carlo 
simulations, where we do know the actual (i.e. input) parameters, do not show correlated 
errors. 

We are measuring the scale height and the velocity dispersions for the tangent point 
gas. The fact that we are sampling a portion of the total gas in the inner galaxy may be 
of concern in that the parameters may not be very representative. The velocity crowding 
at the tangent points means that we are typically looking at ~ 1 kpc of path length. With 
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the present analysis we fit the parameters to intervals of a few hundred parsecs in Galactic 
radius, about 1 kpc in the line of sight distance. The typical number of clouds at tangent 
points is 40 (at longitude / = 35°), so we are not finding the parameters of just a single 
cloud unless the tangent point regions are severely underpopulated. 

4.1. Tangent point analysis 

Since we expect the velocity dispersion to have values of few kms -1 , the tangent point 
emission has to be modeled more carefully than if we were seeking to derive the rotation 
curve. Most of the analyses of the tangent point gas have been done to determine the 
rotation curve of the gas (atomic or molecular). In this paper we are concerned more with 
the determination of velocity dispersion and the scale height of the gas. In principle the 
velocity dispersion and scale height could be estimated simply by looking at the gas at 
forbidden velocities, i.e. V > Vr- This is not very practical however for the following two 
reasons. The first has to do with the calculation of the terminal velocity; most estimates 
involve assumptions about the value of the velocity dispersion, or at least its constancy 
(Gunn et al. 1979), or the size-linewidth relation of the constituent clouds (Clemens 
1985, Liszt et al. 1981). It would be quite circular then to use such a rotation curve to 
calculate velocity dispersions. One could use the rotation curve determined from HI data 
to determine the velocity dispersions for molecular gas. The disadvantage in doing so is 
mainly due to the clumped nature of molecular gas, which means we see only a few clouds 
at forbidden velocities. For this reason the velocity dispersion and scale heights found by 
this method are very sensitive to the terminal velocity. Shane and Bieger-Smith (1966) and 
Burton and Gordon (1978) subtract the equivalent width of the terminal velocity feature, 
which overestimates the velocity dispersion, and underestimates the terminal velocity. 
These methods give values of terminal velocities accurate to a few kms -1 - about the value 
expected for the velocity dispersion. 

4.2. Velocity dispersion 

Various methods have been used to estimate the velocity dispersion of molecular gas. 
These are summarized by Stark & Brand (1989). The values of a v found are also diverse, 
from 4.2 kms -1 to 7.8 kms -1 . The situation is further complicated by the definition 
of velocity dispersion, and its possible confusion with streaming motions. As a working 
definition we assume that any motion that averages out for a distance of AR = 200 — 400 pc 
does not contribute to the velocity dispersion. This is achieved by keeping the parameters 
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(for example rotation velocity 9) stiff over such intervals in radius, while fitting the terminal 
velocity feature (section 2.3). In treatments where the terminal velocity is allowed to vary 
between longitudes quite close to one another ( e.g. figure 3 in Burton & Gordon (1978)) 
smoothing the rotation curve to AR = 200 — 400 pc and calculating the scatter of Vr 
around the smoothed value would be a method equivalent to ours. 

Stark (1984) and Stark and Brand (1989) estimate the velocity dispersion of a 
population of local clouds, at distances of up to 3 kpc from the sun. They find a velocity 
dispersion a v = 7.8^1 kms -1 . Given that they do not distinguish between streaming and 
random motions, and have an operational definition similar to ours; their value of a v found 
in the solar neighborhood is consistent with our value. 

The a v determination likely to be most free of the vagaries of the rotation curve and 
large scale motions is the velocity dispersion derived for high latitude clouds by Blitz et. 
al. (1984) who find a v = 5.7 ± 1.2 km s -1 for a sample of 28 clouds, after excluding a 
'pathological cloud' with velocity — 24.6kms _1 . This is however a measure of the velocity 
dispersion (mostly) in the vertical direction which may be different from the line-of-sight 
velocity dispersion estimated in this paper. Note that these are smaller clouds, and if they 
have a velocity dispersion greater than the galactic plane GMCs (Stark 1983) the vertical 
velocity dispersion value may be still smaller than the velocity dispersion in the plane. 

4.3. Vertical distribution of the gas 

The vertical distribution of the gas is found in this study is consistent with previous 
values. The scale height is found to increase with radius. A linear fit gives a gradient in 
scale height of 5.5pckpc _1 , and extrapolates to give a scale height 58 pc at -R Q . Similar 
values are found by Dame and Thaddeus (1985), Dame et. al. (1987), Sanders et. al. 
(1984), Grabelsky et. al. (1987), and Clemens et. al. (1988). We find from the Monte 
Carlo simulations that the (cloud-cloud) scale height is overestimated if the cloud sizes are 
comparable to the scale height, and because of cloud-cloud shadowing at tangent points. 
The second effect is expected to be smaller at tangent points seen at higher longitude. 
This introduces an error in the gradient of the scale height with radius. The first effect 
(cloud size) will lead to an incorrect constant term in the scale height-radius relation. 
Extrapolating the a z to R we may under or overestimate a z (R ) depending on which of 
these effects dominates. Local determinations of the scale height at Rq for clouds/HII 
regions (Fich & Blitz 1984) give a value of 88 pc. The distinction between the distribution 
of cloud centers and the scale height of the gas may be best made for local clouds. 

The positions of the CO layer centroid, z , are also consistent with the above mentioned 
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studies. The centroid z also shows a dip, going as far as ~ 50 pc (comparable to a z ) south 
of the midplane. Simulations made with the distribution of gas centered at z=0 do not 
show large deviations from z=0 (figure 5, table II). This feature is not seen in the studies 
of the vertical distribution of HI (Celnik et al. 1979). Unfortunately it is difficult to see 
the azimuthal structure of this waviness, since we can only look at the tangent points. 
At any rate this smooth variation argues against massively incorrect velocity-to-distance 
transforms, since that would disturb the smoothness of such a feature if it was there, and 
would require a conspiracy to produce one if it were not. 

4.4. Mid-plane mass density of the disk 

The mid-plane mass density of the disk has been derived making numerous assumptions. 
We assume that the molecular gas layer is in equilibrium, at least on the larger scales. We 
also assume that there is a single population of molecular clouds at least with respect to 
their velocity dispersions. There may be more than one population of molecular clouds 
with different velocity dispersions. The vertical distribution of molecular gas in that case 
will not be a single Gaussian but a sum of more than one Gaussian. One way to test these 
assumptions is to see whether the vertical distribution is indeed Gaussian. We show in a 
subsequent paper that there exists a population of small clouds at high z. The tail of the 
distribution is the inconsistent with a Gaussian distribution. Similar high tails have been 
found for HI (Lockman 1984). Given the 'noise' in the data due to the small number of 
molecular clouds, we do not think a more complicated model can be derived. 

We also assume that the clouds are independent ballistic particles so the number 
of clouds found in any region is a Poissonian variable. Moreover it is assumed that the 
distribution of cloud centers is dictated by the dynamics of the disk alone. In fact the 
occurrence of a cloud in a region of space is not a Poissonian event as the clouds are 
correlated in space and organized into cloud complexes. 

To examine the vertical equilibrium the value of the vertical velocity dispersion should 
be known, something possible only for local molecular clouds. We are always measuring 
the azimuthal velocity dispersion when we look at the tangent points. The gas being a 
dissipational component we expect the velocity dispersion to be isotropic. The assumption 
of the isotropy of the velocity dispersion can be checked by comparing the a v derived for 
high latitude clouds (in part the vertical component), which is 5.7 ± 1.2 km s" 1 (Blitz et 
al. 1984), with the value of the azimuthal component (extrapolated value at the solar 
neighborhood), which we find to be ~ 7.8 km s" 1 . Within the errors this value is consistent 
with the velocity dispersion for the high latitude clouds and also with the expected shape 
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of the velocity ellipsoid. 

The values of the radial (u), azimuthal (v) and the vertical (w) velocity dispersion 
of (late type) stars in the solar neighborhood are related as: u ~ 2w, v ~ \[2w (Mihalas 
& Binney 1981). Given the uncertainty in the estimated velocity dispersion we cannot 
determine if the velocity ellipsoid of the gas has that shape, even locally. Support for the 
isotropy of the velocity dispersions comes from the shape of the velocity ellipsoid of early 
type stars. The O and B stars show fairly isotropic velocity dispersions (Mihalas & Binney 
1981, p 423); (v 2 )/(u 2 ) ~ 1, although the vertical component (w 2 ) is consistently but 
slightly lower. These stars are young enough that they show mainly the kinematics of the 
interstellar medium. 

We would like to check the mass density profile estimated in this study against previous 
measurements; two such check points are the scale length on the exponential disk (of both 
the light and mass distribution), and the mass density in the solar neighborhood. Using 
the values of velocity dispersion and the scale height extrapolated to solar radius we get a 
Po(Ro) = O.34M pc~ 3 ; extrapolation of the exponential profile fitted to individual po(R) 
gives po(-Ro) = 0.2M Q pc -3 with a formal uncertainty of a factor of 2.5. Using the local 
value of vertical velocity dispersion a v (Blitz et al. 1984) gives po(Ro) — 0.18M Q pc~ 3 . 
These results are consistent with the local mass density determinations from stellar 
kinematics p (R ) = 0.2 (Bahcall et al. 1992, Bahcall 1984a, Bahcall 1984b, Bahcall 1984c) 
and Pq(Rq) =0.1 (Bienayme et al. 1987, Kuijken & Gilmore 1989a, Kuijken & Gilmore 
1989b, Kuijken & Gilmore 1989c, Kuijken & Gilmore 1991a, Kuijken & Gilmore 1991b, 
Kuijken 1991). The best fit exponential mass distribution gives a scale length 5kpc, again 
with an uncertainty of a factor of 2. The scale length of the light in the galactic disk has 
been found to be 3.4 kpc (de Vaucouleurs & Pence 1978) and 5 kpc (van der Kruit 1987). 
Both the scale length and the local po(Ro) found in the present study are consistent with 
the values of 5 kpc and 0.1M Q pc~ 3 estimated by Merrifield (1992). 

The systematic errors in this study as we know them are as follows. The vertical 
velocity dispersion may be lower than the azimuthal velocity dispersion used by a factor 
of \/2. This would lower the mass density estimate from the one given here by a factor of 
two. The scale height may be overestimated here because of cloud size and also because 
of shadowing of clouds by one another. This could raise the mass density estimate from 
the one calculated here by a factor of at least 2.8. Also it may lead to an overestimate of 
the scale length of the mass distribution since the shadowing is lower at high longitudes. 
Given the systematic and statistical errors none of the parameters of the mass distribution 
mentioned in the last paragraph is favored over the other 
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5. CONCLUSIONS 

In this work we have modeled the emission from molecular gas at the tangent points in 
the first quadrant of the Galaxy and measured the terminal velocity, the line of sight velocity 
dispersion, the scale height and the deviation from the galactic plane. The modeling takes 
into account emission from a large path along the line of sight to the tangent point. The 
errors in the parameters are estimated by applying the modeling technique to simulations 
of the galactic molecular cloud population. We find 

(1) The terminal velocities are in good agreement with the terminal velocities from 
HI data, indicating that there is molecular gas at most terminal points, and once again 
pointing to the similarity of the kinematics of atomic and molecular gas. 

(2) The velocity dispersion of the molecular gas increases with radius between 
R= 2.5 kpc and 7.5 kpc. A straight line fit extrapolates to give a local value of 7.8 km s -1 . 

(3) The scale height also increases with the galactic radius. Again a linear fit to the 
scale- height as a function of radius gives a local value of 57 pc. 

(4) The centroid of the a dip below the galactic plane with a maximum excursion of 
-51 pc at R/Rq = 0.7. 

(5) The simulations of the distribution of molecular clouds with 'standard' properties 
show morphology similar to the data. Using these simulations to estimate the errors in 
the parameters, we see that the scale height of the gas is measured to be higher than the 
scale height of the cloud centers, depending on the size of the clouds. This effect cannot be 
corrected for without better knowledge of cloud properties. The simulations also show that 
missing gas at the tangent points leads to an overestimate of the velocity dispersion. 

(6) The midplane mass density of the disk shows a decline with the Galactic radius. An 
exponential disk model fitted to the data gives a scale length of 5 kpc. The local midplane 
mass density obtained from extrapolating this model is po(Ro) = 0.2M Q pc~ 3 . Both the 
scale length and the local mass density estimates are uncertain to a factor of two at least. 

It is a pleasure to thank G. R. Knapp for suggesting this problem and for advice 
during all stages of this work. We would like to thank J. Binney, R. Braun, W. B. Burton, 
C. Gammie and D. Spergel for helpful discussions and R. H. Lupton for his very flexible 
graphics package 'SM'. We also thank the refree for comments. This work was supported 
by NSF grant AST89-21700 to Princeton University. 
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6. Appendix A 

This section shows the latitude-velocity maps of l2 CO(J : 1 — > 0) emission at different 
galactic longitudes from the survey of Knapp et al. (1985) along with the best fit models 
for each line of sight. The latitudes have been converted to the height above the midplane 
at tangent points, according to the relation z = b R cos/, and the y axis is labeled as such. 
This conversion holds only for the tangent point emission at extreme positive velocity. 
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Figure 1: The geometry of the tangent point emission. S indicates the position of the sun, 
R is the Galactic radius at the tangent point, and R' is a sub-tangent point. The models of 
the tangent point emission take into account the emission from the annulus between R and 
R' 

Figure 2: Latitude-velocity maps of l2 CO emission at the longitude 1=22 ° (contour levels: 
1.6, 2, 3, 4, 6, 8, 10, 12 K). The best fit model to the tangent point emission is shown 
superposed on the data. The model shows an abrupt cutoff at the velocity VJiaif, because 
the fitting is done only for V > Vhaif; Half being the velocity at which the tangent point 
emission is half its peak value (cf. Section 2.4). 

Figure 3(a): Latitude-velocity maps of 12 CO emission at the longitude 1=39 ° (contour 
levels: 1.6, 2, 3, 4, 6, 8, 10, 12 K) along with the best fit model. The vertical lines show the 
various values of \4aif- The fitting is done for emission at velocities V > Vhaif- 

Figure 3(b): The best fit parameters of the model for the longitude / = 39: the terminal 
velocity Vr, the velocity dispersion a v , the centroid in z - Z and the scale height a z are 
plotted as a function of the cutoff velocity \4aif- Emission at velocities > Vhaif was modeled. 
This figure shows that the parameters of the model are not very sensitive to small changes 
in the boundaries of the region they are fit. 

Figure 4(a), 4(b) and 4(c): Monte-Carlo simulations of b-v emission maps at longitudes 
1=22.02°, 1=35.1° and 1=55.59° respectively. The morphology may be compared to the 
observations at the corresponding longitudes (Figure 12). Ten simulated b-v maps each 
for these and other longitudes (1=25°, 28°, 32°, 39°, 42°, 46°) were modeled to estimate 
the uncertainties in the parameters. The best-fit modela are shown superposed on the 
simulations. The side panels show the z-profile of the tangent point gas, summed over 
V> Vhaif, along with the best fit model. The top panels shows the composite spectrum 
obtained by summing the emission over all observed latitudes along with the model similarly 
summed. 

Figure 5: The average and standard deviations of the best fit parameters for the 
Monte-Carlo simulations of emission seen at longitudes 1=22°, 25°, 28°, 32°, 35°, 39°, 42°, 46° 
and 51°. Ten simulations were done for each longitude. The parameters are - the terminal 
velocity Vr, the velocity dispersion a v , the centroid in z - Z and the scale height a z and 
are plotted against R/Rq = sin(l). The error bars indicate the standard deviation of the 
measured parameters in each sample of ten simulations. Horizontal bars show the input 
parameters. We see that the measured scale height is systematically higher. 

Figure 6: The measured scale height for various simulations for the longitude 1= 35.1° are 
plotted against the input scale height for the distribution of cloud centers (open points). 
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Ten simulations were done for each input scale height, and the error bars indicate the 
variance of the measured scale height in each sample of ten. The contribution due to the 
large sizes of molecular clouds is calculated as <j\ (cloud) = a\ (measured) — (centers). The 
triangles show the values of a 2 z (cloud), which are nearly constant and have an average of 
~ 42 pc (solid line). 

Figure 7: The effect of missing tangent point gas. The lower panel shows the velocity 
crowding dr/dv vs the velocity at the longitude 1=35.1° (dotted line). The dashed 
curve is the convolution of the dr/dv function with a Gaussian of standard deviation 
7, corresponding to a velocity dispersion a v = 7kms~ 1 , and is the expected shape of a 
spectrum. We get the solid curve if there is gas missing from the tangent point. The 
main effect is to broaden the profile. The upper panel shows a similar effect, looking at a 
2-dimensional latitude-velocity map, which is limited in latitude. The thick lines show a 
composite spectrum made by summing the emission in latitude. 

Figure 8: The effect of missing tangent point gas. The best fit parameters for the 
Monte-Carlo simulations of emission seen at longitude 1=35.1° are seen. Triangles represent 
the parameters obtained when gas is missing from the tangent points, open circles represent 
the parameters for the same simulations but with the tangent point gas present. The 
best-fit parameters (especially a v ) show more scatter when there is missing tangent point 
gas. 

Figure 9: The velocity dispersion a z and the terminal velocity Vr for the tangent point gas 
at different Galactic longitudes I (hence Galactic radii R/Rq = sin(Z)). The solid line shows 
the terminal velocity for HI (from Gunn et al. 1979). 

Figure 10: The scale height a z and the position of the centroid of the tangent point gas 
Z at different Galactic longitudes I (hence Galactic radii R/Ro = sin(Z)). The midplane 
deviation from the plane Z at 1=48.6 is comparable to the scale height. 

Figure 11: The midplane mass density po(R) for different Galactic radii. po(R) is calculated 
from the scale height a z and the velocity dispersion a v using equation 10 (section 3.3). The 
solid line shows the best fit exponential disk model with a scale length of 5 kpc. 

Figure 12: Latitude-velocity maps of 12 CO emission at different galactic longitudes (central 
panel). The y-axis shows the height above the plane for the tangent point gas at extreme 
(positive) velocities. Superposed are the best-fit models of the tangent point emission. The 
side panel shows the z-profile of the tangent point gas, summed over V> Vhaif, along with 
the best fit model. The top panel shows the composite spectrum obtained by summing the 
emission over all observed latitudes along with the model similarly summed. The models 
show an abrupt cutoff at velocities Half; because the fitting is done only for V > Vhaif (cf. 
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Section 2.4). The contour levels for the data are: 1.6, 2, 3, 4, 6, 8, 10 and 12 K. The models 
are contoured at 80, 60, 40 and 30 percent of the peak temperature. 
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Table 1. Gas Distribution parameters 
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Table 2. Monte Carlo Simulations 



Longitude Centroid Scale height Terminal Velocity x 2 

Z a z velocity dispersion per degree 

(in ) (pc) (pc) VV-(kms _1 ) a\,(kms _1 ) of freedom 



Input 


0.0 


50.0 




7.0 




22.02 


-4.7 ± 18.0 


71.6 ± 14.5 


112.6 ±5.4 


7.2 ±2.7 
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Table 3. Monte Carlo Simulations for Different Cloud Sizes 
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